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An analysis of glass formation for polymer melts that are diluted by structured molecular additives 
is derived by using the generalized entropy theory, which involves a combination of the Adam- 
Gibbs model and the direct computation of the configurational entropy based on a lattice model 
of polymer melts that includes monomer structural effects. Our computations indicate that the 
plasticization and antiplasticization of polymer melts depend on the molecular properties of the 
additive. Antiplasticization is accompanied by a "toughening" of the glass mixture relative to the 
pure polymer, and this effect is found to occur when the diluents are small species with strongly 
attractive interactions with the polymer matrix. Plasticization leads to a decreased glass transition 
temperature T g and a "softening" of the fragile host polymer in the glass state. Plasticization 
is prompted by small additives with weakly attractive interactions with the polymer matrix. The 
latter situation can lead to phase separation if the attractive interactions are sufficiently strong. The 
shifts in T g of polystyrene diluted by fully flexible short oligomers (up to 20 % mass of diluent) are 
evaluated from the computations, along with the relative changes in the isothermal compressibility 
at T g (a "softening" or "toughing" effect) to characterize the extent to which the additives act 
as antiplasticizers or plasticizers. The theory predicts that a decreased fragility can accompany 
both antiplasticization and plasticization of the glass by molecular additives. The general reduction 
in the T g and fragility of polymers by these molecular additives is rationalized by analyzing the 
influence of the diluent's properties (cohesive energy, chain length, and stiffness) on glass formation in 
diluted polymer melts. The description of glass formation at fixed temperature that is induced upon 
increasing the diluent concentration directly implies the Angell equation (r a ~ Aexp{B/((f>o t p — 
4>p)}) for the structural relaxation time as function of the polymer concentration, and the computed 
"zero mobility concentration" 4>o, p scales linearly with the inverse polymerization index N. 



I. INTRODUCTION 



Recent computations using the generalized entropy 
theory of glass-formation indicate that the strong tem- 
perature dependence of the structural relaxation in poly- 
mer glass-forming liquids is due to packing frustration as- 
sociated with complex monomer shapes, chain connectiv- 
ity, and the stiffness of the side groups and/or backboned 
Weitz and coworkers have likewise established that pack- 
ing efficiency (which is modified by changing rigidity) 
governs fragility in model soft colloidal particle systems, 2 
so that packing efficiency in the glass state seems to be 
generally implicated in affecting fragility. The entropy 
theory of glass formation (ETGF) provides an under- 
standing of many other observed trends in the variation 
of the fragility of polymer glass former with chain stiff- 
ness, cohesive energy, and chain length* 3 - Moreover, the 
calculations reproduce the "internal plasticization"^ 5 - ef- 
fect in poly(ev-olefins) that appears as the length of flexi- 
ble side groups increases^ Our computations also provide 
solid support for the interpretation and empirical gener- 
alization of experimental findings concerning the influ- 
ence of molar mass and pressure on glass formation. The 
theory further implies that profound changes in the prop- 
erties of glass-forming fluids can be induced by modifying 
chain packing through film confinement or the addition of 
diluents, and the theory provides a qualitative framework 
for understanding these changes in terms of variations in 
molecular packing. 



Usually a diluent, i.e., a "plasticizer" , significantly de- 
creases the glass transition temperature T g from that 
of a bulk polymer. Early measurements indicate that 
small molecule diluents (solvents) with a lower T g de- 
press the T g of the host polymer to a greater extent 
than solvents with higher T g &2- Molecular plasticizers 
often impart increased flexibility to the polymer in the 
glass state while also reducing the glass transition tem- 
perature. In contrast to the more widely studied plas- 
ticizers, antiplasticizers also depress T g , but these addi- 
tives increase the stiffness (shear or bulk modulus) of 
polymeric materials in the glass stated— This stiffen- 
ing is of great importance for preservation of biomate- 
rials, such as foods, tissues, and drugs, and for enhanc- 
ing the scratch resistance of polymer films, controlling 
the brittleness and other non-linear mechanical proper- 
ties of polymer materials, etc. Thus, antiplastization 
has many practical applications ! 12 ' 13 Examples of an- 
tiplasticizer/polymer pairs involve tricresyl phosphate in 
polysulfonateii and dibutylphthalate in polycarbonate: 15 

Using a lattice model of polymer glass-formation in 
which an ideal glass transition is identified with the ex- 
trapolated vanishing of the configurational entropy, Di- 
Marzio and coworker s 16 ' 17 conclude that the main effect 
of a solvent additive is to depress the glass transition 
temperature. In particular, they find that the depression 
of T g is a nearly a universal linear function of the plas- 
ticizer mole fraction. Recent molecular dynamics sim- 
ulations have examined the relaxation of bulk polymer- 
solvent mixtures^ and supported polymer films diluted 
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by a good solvent^ near the glass transition temperature. 
The calculations also indicate that T g decreases linearly 
with diluent concentration (volume fraction) over a wide 
range of solvent concentrations (up to 25 These 
earlier theoretical and computational works, however, do 
not explore the influence of diluent structure and thermo- 
dynamic interactions with the polymer, and they also do 
not describe changes in the fragility of glass formation. 

The present work considers oligomeric molecular addi- 
tives and thus is not restricted to "structureless" small 
molecule diluents. We also study systems with a wide 
range of polymer-diluent interaction energies to cap- 
ture general physical trends in polymer/diluent glass- 
formation that are associated with molecular size and 
structure. Our study thus requires the computationally 
non-trivial extension of the generalized entropy theory 
(ETGF) to describe these complex mixtures and thereby 
reveal molecular characteristics that promote plasticiza- 
tion and antiplasticization. The analysis is limited to 
miscible fluid mixtures having a moderate concentration 
of additives (about 10 % mass) where only a single calori- 
metric glass transition exists. The compositional varia- 
tion of T g (and fragility) is evaluated for this relative low 
range of concentrations for miscible binary fluids. 

The generalized entropy theory of Dudowicz, Douglas, 
and Freed (DDF)i i 20 ' 21 focuses mainly on the temper- 
ature range above T g , where a thermodynamic descrip- 
tion of glass-formation is theoretically sensible. The en- 
tropy theory is generalized here to describe glass forma- 
tion in multi-component fluids (e.g., binary mixtures 2 ^), 
provided a single glass transition appears in the multi- 
component fluid. 

The ETGF describes the influence on thermodynamic 
properties of short-range correlations induced by chain 
connectivity, semi-flexibility, and monomer structure. 
The more realistic description of the thermodynamics of 
polymeric systems is achieved using the lattice cluster 
theory (LCT) generalization of the Flory theory for semi- 
flexible chains. This ETGF is combined with Adams- 
Gibbs theor y 22 ' 23 to describe the structural relaxation 
times, 

T Q =r exp{/3A^[ s :/ Sc (T)]}, (1) 

where j3 — 1/ksT, tq is the high temperature limit of 
the relaxation time, A/x is the activation energy at high 
temperatures, and s* is the high temperature limit (ac- 
tually peak value) of the configurational entropy den- 
sity s c (T). Over the temperature range T g < T < Tj, 
the configurational entropy density is found to obey 
the relation s c (T) = sq(1 — Tq/T) that implies, r Q = 
To exp{A/x s*/[sq fcs(T— To)]}, where the fragility param- 
eter at low temperatures is determined from the relation, 

D = [An/(k B ToWJs }. (2) 

DDF use the empirically based relation between A/z and 
Tj, namely A/i/fc^ ~ 6Tj, to compute the fragility pa- 
rameter D without adjustable parameters beyond those 
used in the LCT for the thermodynamics^ (This proce- 



dure greatly simplifies the description of mixtures since 
otherwise the composition dependence of A/i would have 
to be established on a system by system basis.) The 
characteristic temperature Tj is determined as the in- 
flection point in Ts c (T), marking the crossover from the 
range T g < T < Tj, where the product Ts c (T) varies 
linearly with T— ^ to the region of higher tempera- 
tures T > T/ at which the behavior is clearly non-linear. 
Experimental 2 ^— and simulatio n 30 ' 31 data for various 
glass-formers suggest that A/i/fc^ is approximately six 
times the experimental mode coupling temperature T mc , 
which is thus identified in our studies with Tj. The re- 
ported high temperature activation energies A/i for seg- 
mental relaxation of relatively fragile polymers vary from 
12.1 kJ/mol (polybutylene, m = 85) to 17.2 kJ/mol 
(poly(methyl acrylate), m — 102) and tend to be higher 
for strong polymers^ 2 - 

The first calculation of the glass transition tempera- 
ture by DDF 25 uses a Lindemann criterion 3 ^— defined 
in terms of a critical mean square displacement relative 
to the average interparticle separations. However, the re- 
sultant T g are quite close to those obtained from the con- 
ventional (and less controversial) definition of T g as the 
temperature at which the computed structural relaxation 
time equals a value characteristic of glass formation 1 , 
e.g., 100 s. Hence, the present work uses the latter sim- 
pler and common definition. 

Section II focuses on implementing the ETGF to de- 
scribe vitrification, using polystyrene (PS) melts as a typ- 
ical example of a relatively fragile polymer for which the 
influence of different solvents on glass formation is well 
documented. The entropy theory for binary mixtures is 
then applied to examine how the properties, such as the 
size and interaction energies, of diluents added to a fragile 
polymer affect the nature of glass- formation (section III). 
Section IV compares the general trends in glass forma- 
tion for binary fluids and polymer melts. The ETGF is 
then applied in Sect. V to estimate the efficiency param- 
eters for (anti-) plasticizers, where PS is used as the host 
polymer. Section VI describes the compositional varia- 
tion of the structural relaxation times and glass transi- 
tion temperatures of the plasticized polymer melts. The 
temperature and chain length dependence of the extrapo- 
lated "zero mobility" polymer concentration (the analog 
of the Vogel-Fulcher-Tammann (VFT) temperature To 
in glass formation upon cooling) is discussed briefly in 
section VII. 



II. APPLICATION OF ENTROPY THEORY TO 
SPECIES WITH DIFFERENT RIGIDITIES OF 
BACKBONE AND SIDE GROUPS: 
POLYSTYRENE AS A TYPICAL EXAMPLE 

Because cyclic portions of monomer structures have 
not yet been implemented in the entropy theory of glass 
formation, a "compact" structure is introduced to mimic 
the size and shape of a polystyrene (PS) monomer whose 
structure is taken here as comprising two backbone car- 
bons and six side group carbons (see Fig. 1) (each with 
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their bonded hydrogen atoms as united atom groups oc- 
cupying single lattice sites). One bending energy Ej, de- 




FIG. 1: Representation of the styrene monomer as a collection 
of united atom groups (CH n ) used in the lattice cluster theory. 



take the monomer in both the SF polymer and the ad- 
ditive to have the 1-pentene structure (two united atom 
groups in the backbone and three united atom groups 
in the side group per monomer). The computed glass 
transition temperature and fragility of the undiluted SF 
polymer are T ffiP = 321 K and K p = l/D p = 0.30. Al- 
ternatively the fragility of the SF melt is calculated as 
m p = 194 from the relation m = A + In 10 A 2 / D where 
A = log T a (T g ) — log To. For computational purposes, 
To = 10~ 14 s is used as a typical value for the relax- 
ation time limit at high temperatures, and the T g of both 
melts and binary fluids is determined conventionally as 
the temperature for which the computed structural re- 
laxation time is t q = 100 s. 



scribes the local flexibilities of the backbone, and another 
E s describes the side group flexibility. A single effective 
van der Waals energy parameter e is used for all carbon 
atoms to minimize the number of adjustable parameters, 
although the actual energies vary. Hence, the cohesive 
energy e and the cell volume a ce u are the two quantities 
describing intermolecular polymer interactions and de- 
termining the equation of state (EOS) of a polymer melt 
in the entropy theory. The bending energies only exert 
a minor influence on the specific volume, thermal expan- 
sion coefficient, and isothermal compressibility. Thus, 
given the particular molecular topology in Fig. 1, the 
parameters e and a ce u are fit to the specific volumes of 
PS in the range from 350 K to 450 K. Our interest fo- 
cuses on experimental measurements in a AT = 100 K 
region above T g j 38 i 39 The fitting yields the cohesive en- 
ergy e = 275 K and cell volume a ce u = 2.72 A parameters 
for high molar mass PS melts (N = 8000). The backbone 
and side group bending energies of high molar mass PS 
are chosen to reproduce the experimental T g = 375 K 
and fragility m = 13942r— (Experiments have T g vary- 
ing from 360 K to 375 K, and m is reported by different 
authors as m = 124M, 143^, 1464£ The reported 

or extracted VFT parameter D = 1 /K is more scattered, 
as D = 3.83 44 , 5.77i- 8 , and 7.80. 46 ) The determi- 

nation of Eb and E s from experimental data is facilitated 
by the observation that m and T g are linearly correlated 
for a given ratio Eb/E s and also that the side groups in 
PS are stiffer than the backbone. The resultant bending 
energies are determined as E^ = 650 K and E s = 1300 
K, giving the ratio E s /Eb ~ 2. 



III. INFLUENCE OF DILUENT'S PROPERTIES 
ON GLASS FORMATION IN BINARY 
MIXTURES 

The influence of different additives on glass forma- 
tion in fragile high molar mass polymers is examined 
by first considering the SF ("stiff" backbone, "flexible" 
side groups) type polymer, (JV = 8000, E b = 900 K, 
E s = 450 K) diluted by a relatively small amount (10 
% mass) of a fully flexible (Ef, = E s = K) oligomeric 
species with variable properties. The calculations here 



A. Influence of the size of the additive 

The depression in the T g of the mixture relative to the 
undiluted polymer melt is found to increase as the size 
of the additive diminishes and to saturate as the size of 
the diluent becomes large enough. Thus, small oligomeric 
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FIG. 2: Shifts in T g for a SF polymer melt diluted by 10 % 
mass of oligomeric additives as a function their polymerization 
index N. Each curve corresponds to a specific cohesive energy 
e of the diluents. 



species (monomers, dimmers) exert the most pronounced 
depression of T g . Figure 2 depicts the dependence of T g 
of the semi-flexible SF polymer diluted by 10 % mass of 
fully flexible oligomers with variable polymerization in- 
dices N. The fragility of these mixtures also decreases, 
as depicted in Fig. 3. The diminution of the fragility 
K as the additive "shrinks" in size (and N decreases) is 
accompanied by a decrease in the computed fraction of 
excess free volume <p v at T g , a relation previously quan- 
tified by us for polymer melts with variable molar mass 
M. The isothermal compressibility k changes insignifi- 
cantly as N varies, unlike T g (see Fig. 4), and similarly, 
the specific density p of the mixture has no appreciable 
dependence on the sizes of the oligomeric additives (see 
Fig. 5) Thus, decreasing the size of the (fully flexible) 
diluent produces a larger depression of T g and renders 
the mixture a stronger glass-former than the undiluted 
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FIG. 3: Influence of size and cohesive energy of the additives 
on the fragility of a SF polymer melt mixed with 10 % mass 
of the oligomeric diluent. 
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FIG. 4: Ratios of the isothermal compressibility of the mix- 
ture (SF host polymer + 10 % mass of the additive) to the 
pure SF melt evaluated at the T g of the mixture as a function 
of size N of the oligomeric additives. Each curve corresponds 
to the specified cohesive energy e of the diluents. 



polymer. However, varying the diluent's size and flexi- 
bility has only a small influence on the "softness" of the 
material, as measured by its isothermal compressibility k 
and specific density p. The change in the stiffness of the 
antiplasticizer or plasticizer additive is then predicted to 
arise mainly from a change in the average cohesive energy 
density by the additive. 



B. Influence of cohesive energy of the diluent 

The polymer-diluent van der Waals interactions e ps are 
taken as the usual geometric average of the pure compo- 
nent energies, e ps — (eppfiss) 1 / 2 . The computations fix 
the van der Waals interaction strength of the host poly- 
mer at e pp = 200 K 49 , while the diluent energy e ss = e 
is varied. Figures 2 and 3 display the dependence of T g 
and the fragility K, respectively, of mixtures comprising 



FIG. 5: Ratios of specific density of the mixture (SF host 
polymer + 10 % mass of the additive) and the pure SF melt 
at the T g of the mixture as a function of the size N of the 
oligomeric additives. Each curve corresponds to the specific 
cohesive energy e of the diluents. 



the semi-flexible SF polymer diluted by 10 % mass of the 
fully flexible additives with different cohesive energies e. 
Figure 4 presents the analogous dependence of the ratio 
k/k p of the isothermal compressibilities of the mixture 
and the melt evaluated at T g (see the caption). The de- 
pression of T g for the mixture relative to the undiluted 
host polymer diminishes as the cohesive energy of the ad- 
ditive's e increases (at fixed size of the diluent N), while 
the fragility K of the mixture drops relative to the melt. 
Thus, increasing the additive's cohesive energy e (at fixed 
size N) produces a stronger glass-former compared to the 
undiluted melt but diminishes the depression of T g . The 
shift AT g is approximately linear in e at fixed AT, analo- 
gous to the variation of T g with e for polymer melts found 
by us previously. 3 

Recent work by Riggleman et aL 31 indicates that the 
addition of antiplasticizing molecular additives to poly- 
mer materials reduces the compressibility of the mixture 
and correspondingly increases the shear and bulk moduli 
in the glass state, while glass plasticizing additives can 
be defined as molecular additives that have the opposite 
effect. It is important to realize, however, that antiplas- 
ticization normally exerts an opposing influence on the 
shear modulus (high frequency modulus, of course) in 
the melt state, so that the best antiplasticizers of the 
mixture in the glass state are also the best plasticizers 
in the fluid domain.— Thus, discussions of plasticization 
and antiplasticization must also specify the temperature 
regime considered ] 53 ' 54 The present paper strictly refers 
to plasticization and antiplasticization of the glass state, 
and we use the compressibility and density together as 
suitable measures of this effect. The computed compress- 
ibility and density do not exhibit an inversion in relative 
values above and below T g , as observed for the shear 
modulus in some systems? 5 - 3 - - — and thus cannot provide 
reliable information about material "softness" in the melt 
state. Riggleman et al^- have not yet studied the case 
of small molecule additives with weak attractive interac- 
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tions to explore the properties of melts with plasticizing 
additives. 

The shift in the isothermal compressibility (density 
variance) of the diluted polymer depends on the inter- 
molecular interactions in the binary fluid: as the cohe- 
sive energy e = e ss of the additive grows (and, hence, e ps 
for heterocontact interactions also increases), the mix- 
ture's compressibility diminishes. More specifically, our 
model computations indicate that when e = e ss > e pp , 
the compressibility of the mixture diminishes with the 
additive (as does the "softness" of the material), while 
the opposite limit of e = e ss < e pp leads to an increasing 
relative compressibility (the "softness" becomes larger on 
dilution). [The decrease of the compressibility (variance 
of density fluctuations) has been observed in connection 
with the suppression of the /3-relaxation caused by molec- 
ular additives^, an effect suggested to be related to the 
antiplasticization phenomenon; see discussion in ref. by 
Bergquist et aL— ] The influence of e on the relative "soft- 
ness" of the mixtures enables us to classify flexible dilu- 
ents that depress T g into two typical classes: glass plasti- 
cizers and antiplasticizers. If the flexible diluent species 
are relatively small (low N) , the depression of T g is great- 
est (at fixed e). The addition of toluene to PS is proba- 
bly a good example^ T g drops, and the system becomes 
stronger with polymer dilution. When e exceeds e pp of 
the host polymer, the material becomes "toughened" in 
the glass state, a characteristic feature of glass antiplas- 
ticization, while decreasing e leads to "softening" of the 
mixture, a characteristic property of glass plasticization. 
These changes in compressibility as e varies appear for 
all temperatures. The specific density p of the mixture 
(evaluated at T g ) modestly increases as an antiplasticizer 
(large e) is added to the melt, while p tends to diminish 
when a plasticizer (low e) is mixed with the polymer (see 
Fig. 5). This effect, which has often been observed, is 
not by itself sufficient to predict whether an additive is 
an antiplasticizer or not.— ^ Additional information is re- 
quired concerning the density and density variance. As 
e increases (at fixed TV), the fragility K of the mixture 
decreases (Fig. 4) in parallel with the diminishing excess 
free volume fraction cf> v at T g in a similar fashion as found 
previously for polymer melts. 3 

We conclude that the "best" glass plasticizers are small 
flexible species with weak attractive interactions, while 
the "best" antiplasticizers are also compact, but addi- 
tionally have strongly attractive interactions with the 
polymer matrix (polar interactions are often especially 
effective, but the presence of halogen, nitrogen, oxygen 
and sulfur groups can also be effective in creating these 
relatively strong interactions) , as noted from previous ex- 
perimental studies4^£ Increasing e "toughens" the ma- 
terial but diminishes the depression of T g , so additives 
with very attractive interactions are expected to pro- 
duce a small depression in T g (or even an increase of T g ). 
Hence, the size N and the flexibility (as described in the 
next section) of the additives play primary roles in gov- 
erning a mixture's T g , but have little impact on material 
"softness" (as determined by the relative compressibil- 
ity). The cohesive energy e of the diluent changes both 



T g and the "softness" of a material. Dong and Fried 59 
utilize a statistical thermodynamic theory to describe the 
T g s of polymers plasticized by different molecules and ar- 
rive at similar conclusions, i.e., high efficiency plasticizers 
tend to exhibit a weak interactions with the host polymer 
and are small in size relative to the polymer statistical 
segment size. 



C. Influence of semi-flexibility of the additives on 
glass formation 

As the stiffness of the additive molecules increases, the 
depression of T g of the mixture diminishes progressively 
and can even produce an increase in T g . We have com- 
puted the ideal glass transition temperature To, where 
configurational entropy extrapolates to zero, for mixtures 
of semi-flexible SF polymers and 10 % mass of oligomeric 
additives for a range of diluent's sizes (N = 2, 10), co- 
hesive energies (e = 100 K, 400 K), and a wide range of 
bending energies E b (0 K to 2000 K). The ratio of the 
backbone and side group bending energies is kept fixed 
at the same ratio as in the host semi-flexible polymer, 
i.e., Ef,/E s = 2. To varies with Eb in a sigmoidal form, 
leveling off at low and high diluent rigidities (see Fig. 
6). Evidently, the influence of the local stiffness of the 
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FIG. 6: Shifts in the ideal glass transition temperatures To 
for a SF polymer melt diluted by oligomeric additives (10 
% mass) with different sizes N and cohesive energies e as a 
function of the bending energy Eb of the diluent backbone. 
The stiffness of the side groups in the additives is E s = \ j1E\, 
(see text). 



additive on To (and T g ) of the mixture is greater for dilu- 
ents with higher molar masses since the diluents having 
larger monomers contain a greater number of stiff bond 
pairs. (Computations for fully flexible additives display 
a linear correlation between To and T g . The strength of 
this correlation supports the assumption that the trends 
in the variation of To as a function of stiffness also apply 

tO Tg). 
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IV. RELATION BETWEEN THE GENERAL 
TRENDS IN GLASS FORMATION FOR BINARY 
FLUIDS AND POLYMER MELTS: SCALING 
RELATION 

The computed variation of AT g for mixtures with e and 
N (see Fig. 2) can be scaled into the general relation, 

AT g (N,e) = -0 s T g , p ln fe^) , (3) 

where a, b and k are coefficients that are independent of e 
and N. The general scaling relation for AT g of mixtures 
as a function of N and e may be rationalized based on 
the dependence of T g for an (oligomeric) diluent on its 
polymerization index N and the strength of the interac- 
tion potential e. The form of this relation is "derived" 
in section VI which focuses on the concentration depen- 
dence of T g for the mixtures. Equation ([3} successfully 
reproduces the depression of T g when the host high mo- 
lar mass polymer is diluted by different oligomeric fully 
flexible additives. (The caption of Fig. 7 presents the 
numerical parameters providing the best fit of Eq. (J3j) to 
the data obtained using the entropy theory). Equation 
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FIG. 7: General scaling relation for the shift in the glass 
transition temperature AT g for mixtures of a SF polymer melt 
diluted by an oligomeric additives when both the size N and 
interaction energy e of the diluent vary. The relation has 



the form AT a 



-T 9 ,p0sln{(a + b/N)/(l + fee)} with a 



3.11 ± 0.04, b = 4.30 ± 0.09, and fe = (0.53 ± 0.01) x 10" 



([3]) implies that the exponentiated absolute value of the 
relative drop in T g scaled by the factor /(e) = 1 + fee 
should be a linear function of inverse segment size 1/N 
of the additive, where k is a properly chosen constant, 
and all computed data for different N and e from Fig. 2 
can be scaled by an e-dependent linear factor /(e) to fall 
on a single master curve which is a linear function of the 
inverse solvent size 1/N. 

Figure 8 depicts the dependence of the change in 
fragility AK = K — K p of a mixture of the host semi- 
flexible SF polymer with 10 % mass of a fully flexible 
oligomeric species for different diluent sizes and cohe- 
sive energies plotted as a function of the change in ex- 



cess free volume fraction at the respective T g , A<f> v — 
4> v {T g ) — <j>v,p(Tg tP ). Although the data for mixtures 
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Correlation between the change in fragility AK = 
of a SF polymer melt diluted by 10 % mass of an 
oligomeric species and the change in free volume fraction A<f> v 
at the respective T g . 



with different e do not merge onto a single master curve, 
there is a strong correlation between the fragility K and 
the excess free volume fraction <fi v (T g ), a trend qualita- 
tively analogous to the behavior found for the fragility K p 
of polymer melts, establishing a correlation of K p with 
<f> v (Tg iP ) as one of the parameters (N, e, or Eb) varies for 
a given molecular topology. 3 - Thus, the general trends in 
the variation of T g and the fragility for glass formation in 
mixtures are qualitatively similar to the trends revealed 
for pure polymer melts. For instance, there is an evident 
analogy between (a) the influence of the size of diluents 
on T g and the fragility K for the mixtures and on (a') 
the effect of the chain size on T g ^ p and K p for polymer 
melts. Similarly, changes in the cohesive energy affect T g 
and the fragility K identically when (b) e ss of the diluent 
in binary fluids and (b') the cohesive energy e pp in the 
polymer melt increases (or decreases). 



V. PS-DILUENT MIXTURES: THE 
EFFICIENCY OF PLASTICIZATION AND 
ANTIPLASTICIZATION 

Small amounts of a plasticizer lead to a shift in T g of 
the binary fluid that is usually linear in the mass fraction 
4>s of the additive, 



Ao 4> s 



(4) 



where Ao is defined as the "plasticizer efficiency pa- 
rameter" . In particular, Mauritz et al. develop this 
linear relation from diffusion theory and use it to 
model (poly)vinylchoride (PVC) plasticized by di-alkyl 
phthalate o 60 ^ 61 . However, Eq. ((4| fails at higher concen- 
trations (above 40 % mass of diluent) , where the depen- 
dence of Tg on dilution generally becomes nonlinear— 
(see Sect. VI). 
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A series of entropy theory computations for PS (N = 
8000) diluted by % to 20 % mass of an oligomeric 
species with different properties (N, e) reveals whether 
the general trends found for SF polymer systems diluted 
by a fixed amount of additive display analogous behav- 
ior for the plasticizer (or antiplasticizer) efficiency pa- 
rameters Xq. PS is chosen to facilitate some comparison 
with experimental data because PS is one is the most 
extensively studied plasticized polymers. Due to limi- 
tations in modeling the monomer structure of PS (and 
of some diluents as well), the focus in the present work 
is on describing general trends rather than in quantita- 
tively reproducing the experimental data. The plasti- 
cizers are considered as fully flexible oligomeric species, 
i.e., with the styrene subunit structure (Fig. 1). Fig- 
ures 9 to 11 describe the general trends as follows: The 
composition dependence of T g for a mixture of PS and 
% to 20 % mass of diluents is almost linear. In ac- 
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FIG. 9: T a of polystyrene (PS) diluted by either a fully flexible 
pentamer [N = 5) or a dimer (N = 2) with the structure 
of PS but with low and high van der Waals energies e as a 
function of the composition of the mixture. 



cord with results in Sect. Ill, the most efficient plasti- 
cizers (large Ao) are those with the smallest e and small 
sizes. Weak polymer-additive interactions "soften" the 
polymer matrix more than the undiluted material. The 
"compactness" of the additive also is the primary feature 
driving the increase of Ao for antiplasticization. Strongly 
attractive polymer-diluent cross-interactions e ps lead to 
a toughening of the material but diminish the efficiency 
of the depression of T g (see Fig. 9, 11). Thus, the rate of 
decrease of T g with dilution tends to be smaller for an- 
tiplasticized polymers, as demonstrated in Fig. 9. The 
computed Ao vary from 81 to 446, while experiments for 
PS indicate that Ao lies in the range from 241 to 524 de- 
pending on solvent/plasticizeri^. Both plasticizers and 
antiplasticizers reduce the fragility of the host polymer 
(see Fig. 10). The fragility parameter m is also almost a 
linear function of 4> s for relatively low amounts of addi- 
tives (within the diluent range from % to 20 % mass), 
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FIG. 10: Composition dependence ((j> s ) of the fragility K = 
l/D of a PS melt diluted by a fully flexible N=5 and N=2 
oligomeric additive with low and high cohesive energies e. 
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FIG. 11: Composition dependence ((j> s ) of the ratios of 
isothermal compressibility of the mixtures (PS diluted by fully 
flexible N=5- and N=2-mers with low and high cohesive en- 
ergies e) and the undiluted polymer, evaluated at the T g of 
the mixtures. The inset presents the identical relations for 
the specific density ratios at the T g of the mixtures. 



with the slopes 770 varying from 134 to 202 in the model 
computations. 



VI. COMPOSITION PROFILE OF 
CONFIGURATION AL ENTROPY IN MIXTURES: 
T a OF THE PLASTICIZED POLYMER 



We have computed configurational entropies for fully 
miscible binary fluids comprising high molar mass semi- 
flexible polymers and low molar mass diluents over a wide 
range of compositions. Data of this type are available, for 
instance, for PS-plasticizer mixtures where the diluent is 
a solvent, such as toluene, chloroform, carbon disulfide, 
etc. The previous discussion focuses on the range from 
% to 20 % mass of the additives where the composi- 
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tion dependence of T g emerges as linear to good accuracy 
(when the diluents are fully flexible species). To check 
whether the composition profile of T g levels off at high 
diluent content (as observed for many polymer-solvent 
mixtures^ ' 63 i 64 ), we have performed entropy theory com- 
putations for PS diluted by low molar mass (N = 5) 
semi-flexible (Ei, = 600 K, E s — 200 K) oligomers (de- 
noted as species A) to mimic binary fluids in which the 
additive A has a (measurable) glass transition temper- 
ature (T S s = 200 K), which is typically much lower 
than the transition temperature of the host polymer 
{T 9:P — 378 K). The monomers of A have the structure of 
1-pentene and a cohesive energy (e = 200 K) appreciably 
smaller than the one for PS (see section II). Thus, A is 
expected to behave as a typical plasticizer. 

The calculations demonstrate that the inverse of the 
"exponential" factor 1/Q(<f> p ) = (^A^) _1 [s c /s*] in the 
Adam-Gibbs relation is to good accuracy a linear func- 
tion of polymer mass fraction (f> p over the whole compo- 
sition range at fixed T (see Fig. 12). This factor fi com- 




0.2 0.4 0.6 0.8 1 

FIG. 12: Dependence of the inverse "exponential" factor 
1/Q,{4>p) = (/3An)~ 1 [s c /s*] in PS-A mixtures (see text) at 
several fixed temperatures as a function of the mass fraction 
4>p of PS. The inset depicts the temperature dependence of 
the "zero mobility concentration" of the host polymer in PS- 
A mixtures. 



with composition, cj) p as s c ~ a(0o, p — 4> p ). Our compu- 
tations show that the overall composition dependence of 
the inverse "exponential" factor l/tt((f) p ) at fixed temper- 
ature is also linear with slight deviations from linearity 
at high dilutions (see Fig. 12). The linearity of l/fl 
in composition cf> p implies that the computed structural 
relaxation time T a {(j) p ) — tq exp{f2((/> p )} exhibits a com- 
position dependence at fixed temperature that accurately 
follows the Angell equation^ - — 

T a (c/) p ) = T O exp{/3A/i(0 p )[s^((/)p)/s c (0p)]} 

= Aex P {B/(0 O) p-^)}, (7) 

where <p p is the polymer concentration in the mixture, 
0o, p is its "critical" concentration at which the extrapo- 
lated relaxation time diverges (the extrapolated zero mo- 
bility concentration), and A and B are some concentra- 
tion independent constants. 

An analysis of the computed T a (4> p ) over the whole 
composition range yields the variation of T g for PS-A 
mixtures in Fig. 13, where the full composition profile 
of T g is generally nonlinear. The composition dependent 
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FIG. 13: Compositional variation of the glass transition tem- 
perature T g for PS-A mixtures, where 4> p is the mass fraction 
of the polymer. 



bines the concentration dependence of the configurational 
entropy s c of a binary fluid, its high temperature limiting 
value s* , and the high temperature activation energy for 
the PS-A mixtures A/x ~ 6Tj, which is computed using 
the same empirically based relation as for the melts. The 
crossover temperatures Tj as well as the limiting value 
s* are found to be linear functions of the composition 
of the binary mixtures to good accuracy. Therefore, our 
ETGF computations suggest that the high temperature 
activation energy for binary fluids A/i varies linearly with 
additive concentration 

A^(0 S ) = (1 - </> s )A^ p + </> s A^ s , (6) 

For instance, for PS-A mixtures the computed A/x varies 
from 25.6 kJ/mol for PS to 16.5 kJ/mol for A. The con- 
figurational entropy of a binary mixture scales linearly 



efficiency parameter \(4> s ) = —dT g /d<t> s of the plasticizer 
A progressively diminishes as the amount of diluent in- 
creases, a behavior that accords with the observed com- 
positional variation of T g for polystyrene (PS) in 12 dif- 
ferent solvents, where a continuously decreasing negative 
slope dT g /d(j) s is indicated. 6 However, in experiments for 
small amounts of solvent mixed with PS and in our com- 
putations for PS diluted by % to 20 % mass of A, the 
slope X(<j> s ) is constant, so the depression in T g for the 
PS-A mixtures may be evaluated as AT g ~ — A (j) s J& 
Jenckel and Heush have suggested the expressions^ 

T g = (l-(f) s )T g:P + (l) s Tg tS + b(T g:P -Tg tS )(t> s (l-(f> s ), (8) 

for plasticized polymer fluids with b a parameter that 
characterizes the solvent quality of the plasticizer. A fre- 
quently used relation that enables prediction of binary 
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fluid properties from the properties of pure components 
is the Fox equatio n 68 i 69 that is derived by assuming ran- 
dom mixing between the components, equal heat capac- 
ity jumps in the glass transition region, and vanishing 
excess volume of mixing over the entire concentration 
range, to yield 



f 



T, 



g. p 



T 



(9) 



a, * 



Among the other proposed expressions for T g of mixtures 
that also do not require adjustable parameters is the one 
suggested by Fried, 70 



MTg/Tg-p) 



4>sM T g,s/T g , p ) 
{l-<t> s ){T g>a /T g , p ) + 



(10) 



which likewise may be used to evaluate T g of the mixture 
if T 9tP and T ffjS of the pure components are known. A 
variant of this equation 71 contains an extra adjustable 
parameter k\ replacing T g ^ s /T g . p in the denominator of 
Eq. dTOj) to get 



M T g/ T g,p) = 



(j> s \n{T g , s /T g , p ) 
(1 - <t> a )k\ + <t> 3 ' 



(11) 



We use the computed compositional dependence of T g 
for PS-A mixtures and the T g of the pure components 
to test Eqs. ©, (flO"]), and (TTT]). Both equations © and 
(|10p with no adjustable parameters predict a diminishing 
of plasticizer efficiency A as the amount of diluent in the 
mixture increases. The Fried equation does not describe 
our computations well, while the T g prediction by the 
Fox equation closely corresponds to the predictions of 
the entropy theory (see Fig. 13). Our data can be fit 
excellently by the logarithmic relation in Eq. (|lip with 
ki = 0.93 (in contrast, T g ^ s /T gtP ~ 0.53). Analogously, 
the Vogel-Fulcher-Tammann temperature To of the PS-A 
mixtures can be described by the relation, 



ln(T /T , p ) 



<p s ln(T Q . s /T QtP ) 
(1 - 4>s)ki + 4> s 



(12) 



that is identical to Eq. (JTTJ) except for the use of a slightly 
different k\ — 0.89 from the fit (not shown in Fig. 13). 
Since the adjustable parameter k\ is close to unity k\ ~ 1, 
setting ki = 1 in Eq. (JTTJ) produces the simpler relation 
for compositional dependence of T g for the mixtures as, 



rp rp\ — (j) B rp{f) s 

g y,P 3, s ' 



(13) 



and the plasticizer efficiency parameter is 

A(&) = -dTg/d^ = Tl~*' Tf; a HTg, P /T g , ,)• (W) 

When the diluent concentration is small, \(<fi s ) — A(0) = 
Ao, and the glass transition depression AT g = T g —T giP ~ 
— \o4> s can be evaluated using the simpler relation, 



i Tg,p m (^g,p/^3, s)) 



(15) 



which reproduces computations of AT 9 for PS diluted 



by % to 30 % mass of A with relative errors not ex- 
ceeding 5.6 %. The plasticizer efficiency parameter A de- 
rived for A can be generalized because the decrease in T g 
upon dilution occurs for both plasticization and antiplas- 
ticization. Hence, the influence of diluent properties on 
the T g of mixtures may be understood by inserting addi- 
tional relations into Eq. (|15l) describing the dependence 
of the pure component T g a and T g p on various prop- 
erties using relations that have been derived/observed 
when the diluents are oligomeric species. For instance, 
the chain length dependence of T g from experimental 
measurements and from our computations is given by the 
relation l/T g = 1/T g °° + K/N. Additionally, for a given 
molecular topology, T g scales approximately linearly with 
the strength of the van der Waals energy e as T g = A+k e 
(see, for instance, Fig. 3 in Ref. 3). Combining the 
general N- and e-scaling relations for T g with Eq. (fT5")) 
yields Eq. ([3]), which rationalizes the influence of dilu- 
ent properties on T g for binary mixtures and successfully 
reproduces our model computations for a SF polymer di- 
luted by a fixed (small) amount of oligomeric species (see 
section II). Although, these diluents are treated as fully 
flexible and thus have no definite glass transition tem- 
perature, the general trends are excellently described by 
Eq. (fT5|) that has been "derived" from an analysis of the 
compositional variation of T g for PS-A mixtures. 



VII. ZERO MOBILITY CONCENTRATION: 
DEPENDENCE ON TEMPERATURE, 
PRESSURE AND CHAIN LENGTH 

The "zero mobility concentration" (f>o, P generally de- 
pends on temperature in a fashion that can be deduced 
from the entropy theory computations. The inset of 
Fig. 12 depicts the "zero mobility concentration" 0o,p 
for PS-A mixtures at ambient pressure (P ~ 0) as a 
function of temperature for T < Tq. p , where To iP is the 
VFT temperature^ of the undiluted PS polymer melt 
(To p ~ 330 K). The region T > Tq jP is omitted because 
it lacks an extrapolated divergence of the structural re- 
laxation times (no entropy catastrophe) for any compo- 
sition of the PS-A mixture. Thus, <fio, p in the Angell 
equation cannot strictly be interpreted as a "zero mo- 
bility concentration" for T > To iP , as in the case of the 
VFT temperature; To iP only has meaning in the sense of 
extrapolation. 0o. p diminishes monotonically as temper- 
ature decreases, in accord with expectations that if the 
glass formation occurs at fixed temperature by varying 
the composition of binary fluid (e.g., by gradual evapo- 
ration of the solvent) than less diluent <f> s = 1 — <j> p is left 
in the mixture in state Xq when the (extrapolated) di- 
vergence of T a ((p p ) is reached at higher temperature and 
more solvent remains in the fluid in this state Ao at lower 
temperature. 

Our ETGF theory allows determining the pressure de- 
pendence of the "zero mobility concentration". These 
computations are illustrated for PS-A mixtures in the 
temperature range T from 250 K to 310 K and a wide 
pressure range P from MPa to 500 MPa. Basically, </>o, p 
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is linear at low pressures and tend to saturate at high 
pressures. The initial slope d<fio, p /dP for PS-A mixtures 
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FIG. 14: Pressure dependence of "zero mobility concentra- 
tions" 4>o, v °f the polymer in PS-A mixtures at several fixed 
temperatures. Solid lines represent the non-linear fits to Eq. 
(fTo) in the range from MPa to 500 MPa and from 250 K to 
310 K with parameters: 4>o = 1, To — 317.9 K, cpao = 1.335, 
Too = 239.6 K, h = 0.003 MPa" 1 . 

is found to be d(j)o. p /dP ~ —0.001 MPa -1 in the range 
of temperatures from 250 K to 310 K. The temperature 
coefficient d<j) , P /dT ~ 0.004 K -1 from the inset of Fig. 
12, which only slightly varies for different pressures, com- 
bined with the low pressure coefficient d<fio, p /dP from 
Fig. 14, enables us to evaluate the pressure depen- 
dence of the Kauzmann temperature Tr at low pres- 
sures, i.e., the average pressure coefficient of the temper- 
ature at which the configurational entropy extrapolates 
to zero, as dT K /dP = -{d^ p /dP)/{d(j) ,p/dT) ~ 0.25 
K MPa -1 , which is consistent with our direct computa- 
tion of dT K)P /dP ~ 0.30 K MPa" 1 for pure PS at low 
pressures. Moreover, this value is close to the reported 
experimental pressure dependence of the Vogel temper- 
ature dT 0:P /dP ~ 0.32 ± 0.05 K MPa -1 for PS in the 
range from 10 MPa to 80 MPa^. The combined T - P 
dependence of 0o, P for PS-A mixtures (see Fig. 14) can 
be fit well by the expression 



°0,p 



i + (0o - T /T) + hj^ - Too/r)^ 
i + hp 



(16) 



where 4>q, Tq] and 0oo, T^; and h are adjustable pa- 
rameters. The first two parameters are found to equal 
0o — <?Vp,p = 1 an d Tq ~ Tk, p (see also the cap- 
tion to Fig. 14 for numerical values). These parame- 



ters define the temperature dependence of < 



'0,p 



at ambi- 



ent pressure P ~ 0, while the second pair of parameters 
are related to the behavior of 4>o,p m the high pressure 
limit. Lastly, the parameter h -C 1 in Eq. (TT6"|) is as- 
sociated with the low pressure coefficient of the "zero 
mobility concentration", namely d(fio^ p /dP ~ — nh with 

K = 1 + (00 " 0oo) + (Too ~ T )/T. 

The "zero mobility concentration" (j)o, P depends on the 
chain length of the host polymer. The computed 4>o. p at 



fixed temperature for different polymerization indices N 
of the host polymer in PS-A mixtures scales as 4>o,p ~ 
1/N, analogous to the T 9iP ~ 1/A scaling for undiluted 
polymer melts, giving 



°0,p 



w 0,p 



e/N. 



(17) 



The slope 9 of this dependence only changes slightly with 
temperature, as displayed in Fig. 15 where the extrapo- 
lated "zero mobility concentration" of PS is presented 
as a function of the inverse chain length 1/N of the 
host polymer for three different temperatures, T < Xo. p . 
Equation (|15p predicts that the drop in T g upon dilution 
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FIG. 15: Chain length dependence of the "zero mobility con- 
centrations" (f>o, p °f the polymer in PS-A mixtures at several 
fixed temperatures. Solid lines represent the linear regres- 
sions 4>o,p = 0o?p + 0/N with parameters 0g^ p = 0.973 and 
9 = 1.739 (T = 310 K); = 0.902 and 9 = 1.608 (T = 290 
K); </)f£ p = 0.824 and 9 = 1.572 (T = 270 K) 



is proportional to the glass transition temperature T 9yP 
of the host polymer. Hence, the plasticization efficiency 
A(JVj) of a given plasticizer for a lower molar mass host 
polymer should be smaller than X(Nh) for a host polymer 
with higher molar mass by a factor equal to the ratio of 
their glass transition temperatures T giP . To confirm this 
prediction, we have computed AT g of a low molar mass 
PS (N t = 20; T g . p = 361.3 K) diluted by 20 % mass 
of the diluent A which is an oligomer of the high molar 
mass PS {N h = 8000; T 9iP = 378.1 K). The computed 
shift AT g = —44.1 K for the low molar mass PS is close 
to the ATg — —45.6 K predicted using the T Sj p -scaling 
in the relation Eq. (fT5"j) and the computation for the high 
molar mass PS [AT g = —47.7 K). When the molar mass 
of PS drops to N, the predicted shift in T g for the PS-A 
mixture is small because our entropy theory underesti- 
mates the experimentally observed slope of the 1/T g _ p 
vs. 1/N relation for PS systems by a factor of 3. Us- 
ing the experimental dependence of T 9tP on molar mass 
for PS leads to the expectation of a bit larger change 
in the plasticizer efficiency parameter A A/ A ~ 0.1 with 
variation of the molar mass from very high polymeriza- 
tion indices down to Ni — 20 for this example. However, 
the generally favorable qualitative comparison between 
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the observed and calculated AT g for mixtures encourages 
the belief that Eq. (fT5l) is universal enough to capture 
the general trends for glass formation in binary mixtures, 
not only when the properties of diluents vary (as shown 
in section II), but also when some properties the host 
polymer change (e.g., the chain length). 

VIII. CONCLUSIONS 

Systematic computations with the generalized entropy 
theory provide significant insight into the influence of 
additives on glass formation in polymers. Our studies 
reveal general guidelines concerning the molecular basis 
for both plasticization and antiplasticization. The com- 
putations establish the primarily role of diluent's size and 
its interaction energy with the host polymer in creating 
conditions that allow either plasticization or antiplasti- 
cization of the polymeric material. Mixing flexible dilu- 
ent molecules having small size and strong attractive in- 
teractions with the polymer matrix converts the frag- 
ile polymer system into a stronger glass-former with an- 
tiplasticization characteristics, a finding in accord with 
recent simulations of Riggleman et alXS. In contrast, the 
addition small, flexible diluents having weak attractive 
interactions with a host polymer leads to plasticization 
and a reduction in fragility. The computed efficiencies of 
the plasticizer or antiplasticizer are rationalized in terms 
of simple scaling relations that describe the influence 
of solvent size and interactions on the variation of the 
glass transition temperatures with diluent concentration. 
Changes in the isothermal compressibility and specific 
density upon dilution tend to be uniform at all temper- 
atures. For instance, antiplasticized mixtures densify on 
mixing and are less compressible at all temperatures than 
the undiluted polymer melt. However, computations of 
the isothermal compressibility and density rather than 
the bulk modulus cannot reveal the existence of an "an- 
tiplasticization" temperature T a , below which antiplasti- 
cization is observed and above which only plasticization 
appears M 

A simple scaling relation is found to describe the com- 
positional variation of the computed glass transition tem- 
peratures of plasticized polymers in which T g for the 
mixtures is reconstructed from the T g of the pure compo- 
nents. This scaling relation is consistent with experimen- 



tal observations in predicting that the (anti-)plasticizer 
efficiency parameters continuously diminish as the fluid 
becomes more diluted. An analogous relation is deduced 
for the Vogel-Fulcher-Tammann temperatures To as a 
function of composition. The structural relaxation times 
are predicted to obey the Angell equation found pre- 
viously by analyzing experimental data for the viscos- 
ity and electrical conductivity of binary mixtures in the 
supercooled region^ The "zero mobility concentration" 
0o, pi at which the extrapolated relaxation time diverges, 
depends approximately linearly on temperature and on 
the inverse chain length of the host polymer. 

Further research using the entropy theory for mix- 
tures should be aimed at resolving the issue concerning 
whether two glass transition temperatures occur when 
the fluid mixture exhibits phase separation^ In this 
case, the criterion r a ~ 0(100 s) can be applied to deter- 
mine both temperatures simultaneously. (However, the 
criterion cannot be applied if the fluid has two transi- 
tions in the one-phase region). Our computations for 
fully miscible PS-A mixtures do not suggest a sharp vari- 
ation of the configurational entropy of mixtures for tem- 
peratures comparable to the T g of each pure components, 
even when the pure system T g differ by more than 170 
K. 

Another interesting aspect for study is the influ- 
ence of non-Berthelot values of the exchange energy 
( = £pp + £ss — 2 e ps . We have invoked the Bcrthclot 
constraint to restrict the number of interaction energy 
parameters, and lifting this condition should have 
interesting consequences and provide the influence of 
the exchange energy on glass formation in binary fluids 
for a given architecture of the solvent molecule and for 
properties of a host polymer. 
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